The dispersion of growth of matter perturbations in f{R) gravity 
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We study the growth of matter density perturbations 5m for a number of viable f{R) gravity 
models that satisfy both cosmological and local gravity constraints, where the Lagrangian density 
/ is a function of the Ricci scalar R. If the parameter m = Rf,RR/f,R today is larger than the 
order of 10~^, linear perturbations relevant to the matter power spectrum evolve with a growth 
rate s = d In 5™ /din a (a is the scale factor) that is larger than in the ACDM model. We find the 
window in the free parameter space of our models for which spatial dispersion of the growth index 
7o = 7(2 = 0) {z is the redshift) appears in the range of values 0.40 < 70 ;$ 0.55, as well as the 
region in parameter space for which there is essentially no dispersion and 70 converges to values 
around 0.40 < 70 ^ 0.43. These latter values are much lower than in the ACDM model. We show 
that these unusual dispersed or converged spectra are present in most of the viable f{R) models 
with m{z — 0) larger than the order of 10^®. These properties will be essential in the quest for f{R) 
modified gravity models using future high-precision observations and they confirm the possibility to 
distinguish clearly most of these models from the ACDM model. 



I. INTRODUCTION 

The origin of dark energy (DE) responsible for the cos- 
mic acceleration today has been a lasting mystery [l[. 
Although a host of independent observational data have 
supported the existence of DE over the past ten years, no 
strong evidence was found yet implying that dynamical 
DE models are better than a cosmological constant A. A 
first step towards understanding the origin of DE would 
be to detect some clear deviation from the ACDM model 
observationally and experimentally. 

Models such as quintessence [H based on minimally 
coupled scalar fields provide a dynamical equation of 
state of DE different from wde = ~1- Still it is difficult 
to distinguish these models from the ACDM model in 
current observations pertaining to the cosmic expansion 
history only (such as the supernovae la observations). 
Even if we consider the evolution of matter perturba- 
tions S„i in these models, the growth rate of Sm is similar 
to that in the ACDM model. Hence one cannot generally 
expect large differences with the ACDM model at both 
the background and the perturbation levels. 

There is another class of DE models in which gravity 
is modified with respect to General Relativity (GR). The 
simplest one would be the so-called f{R) gravity where 
the Lagrangian density / is a function of the Ricci scalar 
R 0]. The basic idea is that gravity is modified on cos- 
mological scales when R is of the order of {Hq is 
the Hubble parameter today), while Newtonian gravity 
is recovered in the region of high density {R 3> Hq). A 
number of viable /(i?) models have been constructed in 
this spirit 0, i, i, 0, i, 1 EH, [II- Since the law 
of gravity is modified in f{R) models, we can in prin- 
ciple expect large differences with the ACDM model in 
cosmological observations [l^ [ij] and in laboratory tests 
[isl . [l6l . [V^ compared to quintessence models. 



From a cosmological point of view viable f{R) mod- 
els are similar to the ACDM model during the radia- 
tion and deep matter eras ("GR regime"), but impor- 
tant observable deviations from the ACDM model ap- 
pear at late times as the model evolves towards what 
we call here the "scalar-tensor regime" [see below after 
Eq. ([13])]. A useful quantity that characterizes this de- 
viation is m = Rf,RB./f,R [ll, where f^R = df/dR and 
f^B.R = d'^f/dR^. In order to satisfy local gravity con- 
straints we require that m is much smaller than 1 for 
R:$> H^, e.g., m(i?) < 10"^^ for R « 10^ U^^- 
Meanwhile, in order to see appreciable deviation from 
the ACDM model at the background level of cosmolog- 
ical evolution, the parameter m needs to grow to the of 
order at least 0.01-0.1 today. The models proposed in 
Refs. 0, H, H, [13, [l2l constructed to realize this fast 
transition of m. Actually as we will see, the quantity m 
is related to the (critical) scale ~ ~ (3m/_R)^^^ be- 
low which modifications of gravity are felt by the matter 
perturbations. For increasing m and for decreasing i?, 
this critical scale gets larger. 

The modified evolution of the matter density perturba- 
tions 6m provides an important tool to distinguish /(i?) 
models, and generally modified gravity DE models, from 
DE models inside GR and in particular from the ACDM 
model [1^. In fact the effective gravitational "constant" 
Goff which appears in the source term driving the evo- 
lution of matter perturbations can change significantly 
relative to the gravitational constant G in the usual GR 
regime, i.e. Geff — (4/3)G (the "scalar-tensor" regime). 
Then the evolution of perturbations during the matter 
era changes from Sm oc t'^^^ to Sm oc t^V3S-i)/6 ^ where t 
is the cosmic time [1, [l^ . 

A useful way to describe the perturbations is to write 
the growth function s ~ dln^m/dlna as s = (O^)'', 
where flm is the density parameter of non-relativistic 



matter. As well-known one has 70 = 7(2 = 0) ~ 0.55 
p^ . [20j in the ACDM model. It was emphasized that 
while 7 is quasi-constant in standard (non-interacting) 
DE models inside GR with 70 ~ 0.55, this needs not be 
the case in modified gravity models, in particular large 
slopes can appear [2l| (see Refs. [221 for more related 
works) . For the model proposed by Starobinsky Q it was 
found in Ref. [23| that the present value of the growth 
index 70 can be as small as 70 =0.40-0.43 while large 
slopes are obtained. This allows to clearly discriminate 
this model from ACDM. An additional important point 
is whether 70 can exhibit some dispersion (scale depen- 
dence) for viable f{R) DE models. 

The redshift at which the transition of perturbations 
occurs depends on the comoving wavenumber k. It is 
then expected that the resulting matter power spectrum 
has a scale dependence for viable f{R) models. In this 
paper we shall study the dependence of the growth in- 
dex 7 on scales relevant to the linear regime of the matter 
power spectrum. We consider most of viable f{R) models 
proposed in literature to understand general properties of 
the dispersion of perturbations. This analysis will be im- 
portant to distinguish between the f{R) models and the 
ACDM model in future observations of galaxy clustering 
and weak lensing. 



II. COSMOLOGY IN f{R) GRAVITY 
We start with the action 

5 = ^ y d4.TV^/(i?) + Srnig^., , (1) 

where = SttG (G is bare gravitational constant), and 
Sm is a matter action that depends on the metric g^j^ 
and matter fields ^',„. Since we are interested in the 
cosmological evolution at the late epoch, we only consider 
a perfect fluid of non-rclativistic matter with an energy 
density pm- In the following we use the unit = 1, 
but we restore the gravitational constant G when the 
discussion becomes transparent. 

In the flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) spacetimc with scale factor a{t), the variation 
of the action ([T]) leads to the following equations 

3FH^ = p,n + iFR- f )/2-3HF, (2) 
-2FH = pra+F-HF, (3) 

where F = df/dR, H = d/a, and a dot represents a 
derivative with respect to the cosmic time t. The Ricci 
scalar R is expressed by the Hubble parameter H as R = 
6{2H^ + H). In order to study the cosmological dynamics 
in f{R) gravity, it is convenient to introduce the following 
variables 



together with the matter density parameter 

Pn 



3FH^ 



We then obtain the following dynamical equations [J 

— 3X2 + xl — X1X3 , 
X2(2X3 - 4 - Xi) , 

- 2x3{x3 - 2) , 
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where a prime represents a derivative with respect to 
A'^ = Ina, and 

, ^ Rf.RR Rf.R X3 

mir) = — , r = — = — . (9) 

J,R J X2 

One has m = for the ACDM model (/(i?) = R-2A), 
so that the quantity m characterizes the deviation from 
the ACDM model. Since m is a function of r = X3/X2, 
the above dynamical equations are closed. For given 
hmctional forms of f{R), the background cosmological 
dynamics is known by solving Eqs. ([HI)-® with Eq. 
The matter point Pm and the de Sitter point P^s corre- 
spond to [1] 



P„ 



{xi,X2,X3) = 

A -| m(7+10m) 

S im — J- 2(l+m)2 1 



3m 
1+m ' 



l+4m l+4m 
' 2(l-|-m)^ ' 2(l+m) 



WcS 



1+m ■ 



• ^'ds: ixi,X2,X3) = (0,-1,2), 

= 0, Weff = -1 • 

We require that m w to realize the matter era with 
Clm ~ 1 and WcS ~ 0. From the definition of r in 
Eq. ^ we have r = — m — 1 for P„i , so that the mat- 
ter point corresponds to (r, m) ~ (—1,0) in the {r,m) 
plane. Note that the radiation point also exists around 
{r,m) « (—1,0) The de Sitter point P^s correspond 
to the line r = —2. 

Let us next consider linear perturbations about the 
flat FLRW background. In the so-called comoving gauge 
[23 | where the velocity perturbation of non-relativistic 
matter vanishes, the matter perturbation = Sp„i/ Pm 
and the perturbation 6F obey the following equations in 
the Fourier space 



S,n + { 2H + -^j S„i - ^S„ 




-6H'^ + —]SF + 'SHSF + 'SSF 



, (10) 



5F + 3H5F + ( ^ 



l\ 6F 



(11) 



where fc is a comoving wavenumber. We also define 
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which corresponds to the mass squared of the scalar-field 
degree of freedom, the scalaron introduced in [2^ , in the 
region AP ^ R. As we will see later, the quantity m 
remains smaller than the order of 0.1 in the cosmic ex- 
pansion history and the condition ^ i? is largely 
satisfied at high redshifts. 

For cosmologically viable models, the variation of F 
is small (|F| ^ HF) so that the terms including F 
can be neglected in Eqs. pO|) and pT|) . If we neglect 
the oscillating mode of 5F relative to the mode in- 
duced by matter perturbations 5m, it follows that 8F ~ 
Pm<5m/[3(fc^/a2 + from Eq. (HI]). For the modes 

deep inside the Hubble radius {k"^ /a? ^ H^) the domi- 
nant term in the square bracket of Eq. (fTO|) is {k"^ / a?)5F . 
We then obtain the following approximate equation for 
matter perturbations JJ,. ,26i] 



where 



Goff 



G ^k^/a^ + 3M2 



F3(fc2/a2-f M2) 



G 
F 



1 

1+ 



kyja^AP) 
3 l + /c2/(a2M2) 



(13) 

(14) 
(15) 



Here we have restored bare gravitational constant G. An 
equation of a similar type was found for scalar-tensor DE 
models [13] with an essentially massless dilaton field. The 
quantity Goff encodes the modification of gravity in the 
weak-field regime due to the presence of the dilaton in 
the case of scalar-tensor DE models, and of the scalaron 
in our f{R) models. An essential difference is that the 
scalaron mass M in the region of high density can be 
large, inducing a finite range ^ M^^ for the (Yukawa 
type) "fifth-force" . 

Under the linear expansion of the Ricci scalar R in 
the weak gravity background of a spherically symmetric 
spacetime, the effective Newtonian gravitational constant 
is given by [l^, 



G, 



(N) 
eff 



G 
F 



1 



1 



-Mr 



(16) 



where r is the distance from the center of symmetry. 

Poisson's equation in Fourier space with G replaced by 
G^.g'' as given in Eq. ^B)) corresponds to the gravitational 
potential in real space (between two unit masses) V{r) = 
— Gg^^r. It corresponds to the modification of gravity 
in the weak-field regime that is felt by the cosmological 
perturbations. 

For the validity of the linear expansion of R the mass 
M needs to be light such that Mvc <C 1, where rc is 
the radius of a spherically symmetric body. For small 
scales satisfying r <C M^^, the fifth- force reaches its 
maximal value and the modification of gravity reduces 
to a shift of the gravitational constant G — > 4G/(3i^). 
This is the regime we have called here the "scalar-tensor" 
regime. Cosmologically the corresponding shift in the 



scalar-tensor DE models mentioned above depends on 
time and occurs on all scales. Note that in the massive 
regime Mrc ^ 1 the result is no longer valid. It is 
exactly in this regime where the chameleon mechanism 
[28t begins to be at work so that the effective gravita- 
tional constant becomes very close to G to satisfy local 
gravity constraints (as we will see in the next 

section). 

In the region M'^ ^ k'^ /a^ the cosmological effective 
gravitational constant (fTSjl reduces to the form Ggff — 
G/F. Then the evolution of 6^ during the matter era 
is given by Sm oc i^^^ (note that ~ 1 because the 
deviation parameter m is much smaller than 1). In the 
region A/2 < fc2/a2 j^j^yg ^ 4G/(3F), so that the 
matter perturbation evolves as 5m « t(v^-i)/6 during 
the matter era Q. 

The perturbation equation has been derived for 
sub-horizon modes under the neglect of the oscillating 
mode (such as the term 5F). For cosmologically viable 
models, the solutions obtained by solving Eg. (|13D agrees 
well with full numerical solutions [l^, |29| . l30l|. In our 
numerical simulations, however, we shall solve the full 
perturbation equations ((TO)) and pT|) together with the 
background equations ([S])-®, without relying on the ap- 
proximate equation p3p . For the numerical integration 
it is convenient to rewrite Eqs. PT!]) and pT|) in the fol- 
lowing forms (lo| 



5,„ 



2^3 - ^a;i ) 5'm - X-i)5r, 



fc2 



- 6 + 3a:? - 3a:'i - 3a;i(a;3 - \)\ 5F 



-h3(-2a;i + .T3 - \)5F' + 2,5F" 



5F" + (1 - 2x1 + xs)5F' 



2X3 



2X3 



Xl{x-i + 1) — x'l + x\ 



= (1 - Xi - X2 - X3)5m - Xi5m 



(17) 

5F 
(18) 



where 5F = 5F/F, and the new variable .T4 = aH satis- 
fies 

4 = (X3 - 1)0:4 . (19) 
The growth index 7 of matter perturbations is defined 

as 



s = {nm)\ (20) 
where s = dln(5,„/dlna = 5'm/ 5m and f2„i is given by 

Pn 



3iJ2 



Fn„ 



(21) 



This choice of flm comes from rewriting Eq. ([2]) in the 
form 3i?2 = Pm+PDE, where poE = {FR- f)/2-3HF+ 
3i/2(i _ F) [i,!!!]. For viable f{R) models the quantity 
F approaches 1 in the asymptotic past because they are 
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similar to the ACDM model (as we will see in the next 
section). Defining as well as flm in the above way, 
the Friedmann equations recast in their usual General 
Relativistic form for dust-like matter and DE. 



III. VIABLE f{R) MODELS 

Let us briefly review viable f{R) models that can sat- 
isfy both cosmological and local gravity constraints. We 
focus on the models in which cosmological solutions have 
a late-time de Sitter attractor at i? = i?i (> 0). For the 
viability of f{R) models the following conditions need to 
be satisfied. 

• (i) Jm > for i? > i?i (> 0). This is required to 
avoid anti-gravity. 

• (ii) f.RR > for i? > This is required for the 
stability of cosmological perturbations [l^ , for the 
presence of a matter era |3ll|. and for consistency 
with local gravity tests [1^ • 

• (ii) f{R) ^ i? - 2A for i? > i?o, where Ro is the 
Ricci scalar today. This is required for consistency 
with local gravity tests lUland for the presence of 
radiation and matter eras J. 



• (iv) < m{r ~ —2) < 1. This is required for the 
stability of the late-time de Sitter point [3, [13, HI] . 

The conditions (i) and (ii) mean that m = Rf.RR/f,R > 
for R > Ri. The trajectories starting from {r,m) sa 
(— 1,+0) to the de Sitter point on the line r = —2, < 
m < 1 are acceptable. In other words, the deviation 
parameter m is initially small (0 < m ^ 1) so that the 
model is close to the ACDM model during the radiation 
and deep matter eras. The deviation from the ACDM 
model can be important at the late epoch. Depending on 
the models of f{R) gravity, the parameter m can grow 
as large as to the order of 0.1 today. 

We also require the condition < m ^ 1 in the region 
R^ Ro for consistency with local gravity constraints. In 
this case the mass squared in Eq. ([12]) becomes large 
in the region of high density to avoid the propagation 
of the fifth force. It is then possible for f{R) models 
to satisfy local gravity constraints under the chameleon 
mechanism [2^ . We introduce a new metric variable g^i^ 
and a scalar field 0, as = ^^'^ ^ ~ ^3/2 In ip, 

where tp = F (R) . The action in the Einstein frame is 
then given by [3J| 



S ^ J d-^xy^ \^R/2 - {V(j)f/2 ~ 1/(0) 
+^™(gM-e2W,*„0, 



where 



_ RW^ - f 
V6 2^-2 



(22) 



(23) 



The scalar field degree of freedom has a constant 
coupling Q with non-rclativistic matter in the Einstein 
frame. 

In a spherically symmetric spacetimc of the Minkowski 
background the field obeys the following equation in the 
Einstein frame 



d^ 
df2 



2d0 
f dr 



dFeff 



(24) 



where f is the distance from the center of symmetry and 

V,s{cj,)^V{cj,)+eQ'^p*. (25) 

Here p* = e^'^'^Pm is a conserved matter density in the 
Einstein frame [11]. For f{R) models {Q = —l/\/6) 
the effective potential 14^(0) possesses a minimum for 
Vtj,{4') > 0. For a spherically symmetric body with con- 
stant densities pA and pB inside and outside the star, 
the effective potential has two minima at the field val- 
ues 0A and 0B satisfying the conditions VcS,4>i4'A) = 
and VcS,^{(I)b) ~ 0, respectively, with mass squared 
m\ = V"cff,00(0A) and m% = VcS. d,d,(<t'B)- We define 
the so-called thin shell parameter [28j 



eth 



(26) 



where $c is the gravitational potential at the surface of 
the body (f = fc). If the field is sufficiently heavy such 
that mATc 3> 1 and if the body has a thin-shell in the 
region ri < r < rc with Afc = fc — fi <C ^c, the thin- 
sh ell p arameter is given by eth — ^^cl^c + XjimATc) ^ 
1 [33 • The effective coupling between non- relativistic 
matter and the field is Qcff — 3(5eth, whose strength 
can be much smaller than 1 for eth 1- 

The tightest bound on eth comes from the solar sys- 
tem test of the violation of equivalence principle for the 
accelerations of Earth and Moon toward Sun [2I] . This 
constraint is given by [3^ 



eth,© < 



X 10-V]Q], 



(27) 



where eth,© is the thin-shell parameter for Earth. Since 
the gravitational potential for Earth is ^c,© = 7.0 x 
10~^", the condition ([TT)) translates into 



bB,®\ < 3.7 X 10- 



(28) 



where we have used ]0s,©] ^ ]</'A,ffl]- For cosmologically 
viable f(R) models, ]0b,©] is roughly the same order as 
the deviation parameter m(i?B) at the Ricci scalar Rb 
@ (as we will see below). Hence the condition m{RB) ^ 
10~^^ needs to be satisfied in the region around Earth 
(in which Rb ^ Rq). Cosmologically this means that 
the parameter m is very much smaller than 1 in radiation 
and deep matter eras. 

We consider f{R) models which can be consistent with 
both cosmological and local gravity constraints. They 
can all be written in the form: 



fiR)^R-XRcfiix), 



X = R/Rc , (29) 
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where Rc (> 0) defines a characteristic value of the Ricci 
scalar R and A is some positive free parameter. 
Wc will study the following models: 

• (A) (0<p<l), 

• (B) /i(a;) = + 1) (n > 0) , 
. (C) A (x) = !-(! + {n > 0) , 
. (D)/l(a;) = l-e-^ 

• (E) /i(.t) = tanh(a;) . 

For A of the order of unity, Rc roughly corresponds to 
the scale of the cosmological Ricci scalar Rq today. 

The model (A) is characterized by m = p{r + l)/r, 
which behaves as m ~ p{—r — 1) during radiation and 
deep matter eras (r ~ — 1) [1]. In the regime R ^ Rc 
we have ~ ~{V6/2)m/{l — p) with m ~ \p{l — 
p){R/ RcY^^ ■ In order to satisfy the condition the 
|)^ameter p is constrained to be very small: p < 3x 10~^° 

In the regime R ^ Rc the models (B) and (C), pro- 
posed by Hu and Sawicki and Starobinsky [H respec- 
tively, behave as 



f{R)^R~XRc [1-iR/Rc)-^"] 



(30) 



which corresponds to 
m(r) = C(-r-l)2"+i, C = 2n(2n + 1)/A^" , (31) 

with the field value (f> ~ ~{\/6/2)m,/{2n+ 1). Because of 
the presence of the power 2n-|-l in Eq. pTjl . the condition 
(PS)) can be satisfied even for n and A of the order of 
unity. In fact the models (B) and (G) are consistent with 
the constraint for n > 0.9 [131 • these models 
the deviation parameter m can grow to the order of 0.1 
today. 

The models (D) and (E), proposed by Linder [l3| and 
Tsujikawa respectively, have vanishingly small m in 
the region R ^ Rc, but m can grow to C(O.l) once R 
decreases to the order of Rc (see Ref. Q for a similar 
model). In the model (D) one has m ~ A(i?/i?c)e^^/^'= 
for R > Rc, where R/Rc ~ log{- y^c/j / \) . The field 
value 4>B can be derived by solving T4ff,0(0B) = with 
p* ~ p B , wh ich gives i? ~ pB- We then find that 
(t)B ^ -^/3/2Xe-P"^^"'. Since \Rc is of the order of the 

cosmological density p^^"^ ~ 10~^^g/cm'^ today, we have 
(j)B ~ — Ae"^'' ^ by using the dark matter/baryonic den- 
sity pb — 10^^^ g/cm^ in our galaxy. As we will see later 
the existence of a stable de Sitter point demands A > 1, 
under which the constraint ((28)) is well satisfied. The 
model (E) also has a similar property. Thus the models 
(D) and (E) are consistent with local gravity constraints 
for A required for viable cosmology. 



IV. GROWTH INDICES OF MATTER 
PERTURBATIONS 

In this section we study the growth of matter pertur- 
bations for the viable f{R) models (A)-(E). We are inter- 
ested in the wavenumbers k relevant to the galaxy power 



spectrum [37 1 



0.01 /iMpc"^ < fc < 0.2 /iMpc" 



(32) 



where h = 0.72 ± 0.08 corresponds to the uncertainty of 
the Hubble parameter today [s^. The scales ([32)1 are in 
the linear regime of perturbations. For k = 0.2 h Mpc~^, 
non-linear effects are still small so that results in the 
linear regime can be related to observations. Non-linear 
effects increase as we go to smaller scales. On the other 
hand we should remember that observations on the large 
scale around k ^ 0.01 hMpc~^ are not so accurate but 
will be improved in the future. 

We recall that the transition from the "GR regime" to 
the "scalar-tensor regime" occurs at = k^/a^. Using 
Eq this translates into 



m « {aH/kf 



(33) 



This expresses in terms of the quantity m the fact that 
cosmic scales smaller than ^ M^^ will be affected by 
modifications of gravity. For viable f{R) models we have 
presented in the previous section, the deviation param- 
eter m increases from the matter era to the accelerated 
epoch. For larger k (i.e. for smaller scales) the transi- 
tion occurs earlier. The upper bound of the wavenum- 
ber in Eq. ([5^ corresponds to fc ~ 600aoi/oi where the 
subscript "0" represents present quantities. We are in- 
terested in the case where the transition to the scalar- 
tensor regime occurred by the present epoch (the redshift 
2 = 0). This then gives the following condition 



m{z = 0) > 3 X 10" 



(34) 



If m{z — 0) < 3 X 10"*^ then the linear perturbations 
have been always in the GR regime in the past, so that 
the models are not distinguished from the ACDM model. 
We caution that the bound (|34p is relaxed for non-linear 
perturbations with k > 0.2 /iMpc^^, but the linear anal- 
ysis is not valid in such cases. 

In our numerical simulations we identify the present 



epoch to be 



(0) 



0.28. 



A. Model (A) 



Let us first consider the model (A). In this model the 
deviation parameter m = p{r + l)/r corresponds to m 
p/2 at the de Sitter point (r = —2), which means that 
m{z = 0) is of the order of p. Hence the condition ((M]) 
for the occurrence of the transition to the scalar-tensor 
regime corresponds to 



p > 6 X 10" 



(35) 
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Figure 1: The growth indices 70 today versus the parameter 
p in the model (A) for four different values of k. Under the 
local gravity bound p < 3 x 10"^", the deviation of 70 from 
that in the ACDM model (70 — 0.55) cannot be seen for these 
wavenumbers. 



In Fig. [T] we plot the growth indices 70 of matter per- 
turbations today for four different wavenumbers k. When 
p > 6 X 10^^ the deviation from the value 70 ~ 0.55 
of the ACDM model can be clearly seen for the modes 
k > 0.2 h Mpc^^. The dispersion of 70 with respect to 
the wavenumbers k is especially significant for 10^^ < 
p < 10"^, whereas 70 converges to a value around 0.4 for 
P > 10"^ 

We recall that local gravity constraints give the bound 
p < 3 X 10~^°, which is not compatible with the condi- 
tion (|35p . Under this bound the growth indices are very 
close to the ACDM value 70 ~ 0.55 for the wavenumbers 
Hence the model (A) cannot be distinguished from 
the ACDM model as long as local gravity constraints are 
respected. 



B. Models (B) and (C) 

We shall proceed to the models (B) and (C). In the 
region R ^ Rc these models can be described by the 
TO(r) curve given in Eq. (|31[) . In the deep matter era 
(r « —1) the deviation parameter m gets smaller for 
increasing n because of the larger power-law index 2n + l 
in Eq. (|3ip . For increasing A we also have smaller m. 

In the model (B) the stability of the late-time de Sitter 
point requires that p6j 

2xrf"-(2n-l)(2n-|-4)x^" + (27i-l)(2n-2) > 0, (36) 

where Xd = Ri/Rc- 



The parameter 



^2 



2xf 



2n) 



(37) 



has a lower bound determined by the condition (j36p . 
When n = 1, for example, one has > \/3 and 
A > 8^3/9. 

Similarly the model (C) satisfies 

(1 + a;^)"+2 > 1 + (n + 2)xl + {n + l)(27i + l)x\, (38) 



with 



2[(1 



-^-l-{n+l)xl] 



(39) 



When 71 = 1 we have Xd > \/3 and A > 8V3/9, which 
is the same as in the model (B). For general n, however, 
the bounds on A in the model (C) are not identical to 
those in the model (B). The minimum values of A are of 
the order of unity in both models. 

At the de Sitter point the model ([3T|l gives m{r ~ 
-2) = C = 2n(2n + 1)/A^", so that m(z = 0) can be as 
large as 0(1) for n, A of the order of unity. Numerically 
we find that the deviation parameter m{z = 0) in the 
models (B) and (C) is typically smaller than that in the 
model ([5T1) . but still m{z = 0) can be of the order of 
0.1. The deviation parameter m needs to be very much 
smaller than 1 in the region of high density (i? ^ Re) for 
consistency with local gravity constraints. 

If the transition characterized by the condition ([55]) 
occurs during the deep matter era (2 ^ 1), one can 
estimate the critical redshift Zc at the transition point. 
We use the asymptotic forms m ~ C{~r — l)2"+i and 
r ~ — 1 — XRc/R as well as the approximate relations 
~ H^nin\l + zf and R ~ 3H^. The present value 
of pde may be approximated as p^g w \Rc/2. Hence we 

have that Ai?c ~ SiJgfipEi where ri^E is the DE density 
parameter today. Then the condition (j33p translates into 
the critical redshift 



QoHq 



2n(2n-f 1) (2f}[5'^)2"+i 



X2n 



(0) 



2(n+l) 



-1. (40) 



For n = 1, A = 3, fc = 300aoffo, and fJ^^^ = 0.28 in the 
model (C) the numerical value for the critical redshift is 
Zc = 4.5, which shows good agreement with the analytical 
value estimated by Eq. ([40]) . We caution, however, that 
Eq. (PD|) begins to lose its accuracy for Zc close to 1. 

We recall that local gravity constraints give the bound 
n > 0.9 for both the models (B) and (C). Meanwhile the 
conditions ((HI)]) - ([55)1 provide lower bounds on A for each 
71 (A > 1.54 for 71 = 1 in both models). In Fig. [2] we 
plot the evolution of the growth indices 7 in the model 
(B) with 71 = 1 and A — 1.55 for a number of different 
wavenumbers. We find a degeneracy of the present value 
of 7 around 70 — 0.41 independent of the scales of our 
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Figure 2: The evolution of 7 versus the redshift 2 in the model 
(B) with n = 1 and A — 1.55 for four different values of k. In 
this case the dispersion of 7 with respect to k is very small. It 
is nearly absent for scales k > 0.033 h Mpc~^, so these scales 
have reached today the asymptotic regime k S> aM. 



interest. In this case the transition redshift corresponds 
to Zc = 5.2 and Zc = 2.7 for the modes k = 0.1 hMpc~^ 
and k = 0.01 hMpc~^ , respectively. At the present epoch 
these modes are in the "scalar-tensor" regime with simi- 
lar growth indices. 

Equation (|40p shows that Zc gets smaller for increas- 
ing n. In Fig. [3] we show the present values of 7 ver- 
sus n in the model (B) with A = 1.55 for four differ- 
ent wavenumbers. We find that 70 has a scale depen- 
dence in the region 0.42 < 70 ^ 0.55 for 2 < n < 7, 
while 7o is degenerate around 0.41 for n close to 1. This 
reflects the fact that, for larger n, the transition red- 
shift Zc gets smaller. The growth indices are strongly 
dispersed if the mode k = 0.2/iMpc~^ crossed the tran- 
sition point at Zc > 0{1) and the mode k = 0.01 /iMpc~^ 
has marginally entered (or has not entered) the scalar- 
tensor regime by today. Since Zc decreases for increasing 
A from Eq. it is expected that the scale dependence 
of 7o can appear for larger A than in the case shown in 
Fig. [3] (for fixed n) . In fact this behavior is clearly seen 
in the numerical simulation of Fig. 31 which shows that 
in the model (B) with n = 1 the dispersion of 70 occurs 
for 0.5 < log;^Q A < 2.5. If A > 10^, 70 converges to the 
ACDM value ~ 0.55 because the modes ([5^ have not 
entered the scalar-tensor regime by today. 

We have also carried out numerical simulations for the 
model (C) and found that the evolution of 7 is very simi- 
lar to that in the model (B) for the same values of n and 
A. Let us consider the parameter regions of {n, A) for the 
models (B) and (C) in which the dispersion of 70 occurs 
for the wavenumbers ([5^ . We can divide the (n, A) plane 
in three regions: 



Figure 3: The growth indices 70 today versus n in the model 
(B) with A = 1.55 for four different values of k. The dispersion 
of 70 occurs in the region 0.42 < 70 0.55 for 2 < n < 7. 
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Figure 4: The growth indices 70 today versus A in the model 
(B) with n = \ for four different values of k. The dispersion 
of 70 appears for 0.5 < logj^Q A < 2.5. 



• (i) All modes have the values of 70 close to the 
ACDM value: 70 = 0.55, i.e. 0.53 < 70 < 0.55. 

• (ii) All modes have the values of 70 close to the 
value in the range 0.40 < 70 0.43. 

• (iii) The values of 70 are dispersed in the range 
0-40 < 70 < 0.55. 

We recall that the first case arises when all scales under 
consideration are close to the asymptotic regime for scales 
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larger today than the range of the "fifth-force". The 
second case corresponds to the opposite situation. In the 
third case some of the scales belong to the intermediate 
regime [39| . To find out accurately when the asymptotic 
regimes are reached, and what are the values of 70 in 
the intermediate regime, one has to resort to numerical 
calculations. 

The region (i) is characterized by the opposite of the 
inequality ((341), i.e. m{z = 0) < 3 x 10"''. This corre- 
sponds to the case in which n and A take large values so 
that m is suppressed. The border between (i) and (iii) 
is determined by the condition m(z = 0) « 3 x 10~^. 
The region (ii) corresponds to small values of n and A, 
as in the numerical simulation of Fig. [21 In this case the 
mode k = 0.01 hMpc~^ at least entered the scalar-tensor 
regime for Zc > 0(1)- 

The regions (i), (ii). (iii) can be found by solving per- 
turbation equations numerically. Note that we also have 
the local gravity constraint n > 0.9 as well as the condi- 
tions (IMl) and dMl) with §7^ and ^ coming from the 
stability of the late-time de Sitter point. In Fig. [5] wc 
illustrate the regions (i), (ii), (iii) for the models (B) and 
(C), which are quite similar in both models. The parame- 
ter space for n < 3 and A = 0(1) is dominated by cither 
the region (ii) or the region (iii). These unusual con- 
verged or dispersed spectra can be useful to distinguish 
the f{R) gravity from the ACDM model. 



C. Models (D) and (E) 



by 



The deviation parameter m in the model (D) is given 



Axe ^ 
1 - Ae-=^ 



(41) 



In the region R ^ Rc we have that m ~ Axe~^, which 
means that m rapidly decreases as we go back to the 
past. In the asymptotic past the model (E) has a similar 
dependence m ~ 8Xxe~^^. In both models the parameter 
r behaves as r ~ — 1 — X/x for R » Re- 

For the model (D) the Ricci scalar at the de Sitter 
point {xd = Ri/Rc) is determined by A, as 



A 



2 - {2 + Xd)e~^'^ 



(42) 



From Eqs. (^1]) and ((^ we find that the stabihty con- 
dition m{Ri) < 1 is satisfied for Xd > 0. It then follows 
from Eq. ([1^ that A is bounded to be 



A > 1 



(43) 



If the crossing AP — k'^/a^ occurs during the matter 
era, the transition redshift Zc for the model (D) can be 
estimated as 



2a 



(0) 
DE 



A2(l+Ze)- 



■ exp 



n'^hii + z.f 



2n 



(0) 
DE 



/ k 



(44) 




Figure 5: The regions (i), (ii) and (iii) for the model (B) 
(top) and for the model (C) (bottom). The three regions in 
the model (B) are similar to those in the model (C). 



The redshift Zc gets larger for increasing k and for de- 
creasing A. If fc = SOOaoHo and A = 1.1 we have Zc = 3.0 
from the estimation (|44|l . This is slightly different from 
the numerical value Zc — 2.7 because the transition point 
is close to the onset of the cosmic acceleration. 

In Fig. [6] we plot the growth indices 70 today versus 
A for four different wavenumbers. If A is close to 1 then 
0.40 < 7o < 0.42, so that the dispersion of 70 is weak. 
The dispersion begins to appear for A > 2. This is associ- 
ated with the fact that the transition redshift gets smaller 
for increasing A. If the condition m{z = 0) ^ 3 x 10^^ 
is satisfied, the transition does not occur by today so 
that 7o is close to the ACDM value 0.55 for the modes 
(15^ . Numerically the present value of x is found to 
be xq w 2.2A. Plugging this into Eq. ([1T|) . we find 
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The dc Sitter point is determined by the relation 



A = 



2sinh(xc;) cosh(a::d) — Xd 



(46) 



From the stabihty of the de Sitter point we require that 

ii 



A > 0.905, Xd> 0.920. 



(47) 



As in the model (D) the numerical value of Xq is about 
xq w 2.2A. It then follows from Eq. pS)) that the condi- 
tion m{z = 0) < 3 X 10^^ corresponds to A > 4, in which 
case 7o is degenerate to 70 — 0.55 for the modes ((32l) . In 
Fig.[7]the dispersion of 70 can be seen for A < 4. When A 
is close to the minimum value 0.905, the growth indices 
are almost degenerate in the range 0.42 < 70 ^ 0.43. 



Figure 6: The growth indices 70 today versus A in the model 
(D) for four different values of k. We note that the transition 
is much faster in terms of A than in model (B) shown in Fig.|4] 
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Figure 7: The growth indices 70 today versus A in the model 
(E) for four different values of fc. It is seen that the transition 
for this model is slightly more rapid in terms of A than in 
model (D) shown in Fig. El 



that the condition ra{z = 0) < 3 x 10^^ translates into 
A > 8. This shows that the dispersion of 70 in the range 
0-40 < 70 < 0.55 occurs for 2 < A < 8. This can be con- 
firmed in the numerical simulation of Fig. [HI For A > 8, 
7o converges to the ACDM value ~ 0.55. 
For the model (E) we have 



2Aa;tanh(a;)[l - tanh^(.T)] 
1 - A[l - tanh^(a;)] 



(45) 



V. CONCLUSIONS 

In this paper we have studied the dispersion of the 
growth index 7 of matter perturbations in /(i?) gravity 
models. We focused on a number of viable /(i?) dark 
energy models proposed in the literature that can satisfy 
cosmological and local gravity constraints. While these 
models are close to the ACDM model in the asymptotic 
past, a deviation from the ACDM model appears at late 
times. A useful quantity that characterizes this deviation 
is given by to = Rf.RRjf,R- This quantity needs to be 
very much smaller than 1 during the deep matter era for 
consistency with local gravity constraints, but a growth 
of TO to the present value of the order of 0. 1 can be allowed 
depending on the f{R) models. 

The transition of matter perturbations from the GR 
regime to the scalar-tensor regime occurs at the epoch 
characterized by the condition to w (aH/k)'^. For the 
wavenumbers k relevant for the observable range of the 
linear matter power spectrum (0.01/iMpc~^ ^ k < 
0.2/iMpc^^), we require that m{z = 0) > 3 x 10"^ for 
the occurrence of such a transition by today. For the 
model (A) this requirement is not compatible with lo- 
cal gravity constraints and hence this model cannot be 
distinguished from the ACDM model. 

The models (B) and (C) allow for a rapid growth of 
TO from the region i? > iJg [with m{R) < 10~^^] to the 
region R c::^ H§ [with m{R) = 0(0.1)]. When n < 3 
and A = 0(1) we find two distinct regions widely spread 
in the (n. A) plane: the region (ii) in which the present 
growth indices 70 almost converge to the values around 
0-40 ^ 7o ^ 0.43 and the region (iii) in which 70 are 
dispersed around 0.40 < 70 < 0.55. In the first region 
there is essentially no spatial dispersion of 70, in contrast 
to the second region. These results are summarized in 
Fig.EJfor the models (B) and (C). 

The models (D) and (E) give rise to an even faster 
evolution of to compared to the models (B) and (C). The 
evolution of 7 depends on the single parameter A. For 
the models (D) and (E) the dispersion of 70 in the region 
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0-40 < 70 < 0.55 is found for 2 < A < 8 and 1.5 < A < 
4, respectively. The growth indices converge to values 
around 0.40 < 7o < 0.43 for 1 < A < 2 [model (D)] and 
for 0.905 < A < 1.5 [model (E)]. 

We have thus shown that the dispersed or converged 
growth indices with 70 smaller than 0.55 are present in 
viable f{R) models with m(z = 0) > 3 x 10~^. If future 
observations detect such unusually small values of 70 , this 
can be a smoking gun for f{R) models. The presence 
of some dispersion in future observations could be an 
additional evidence for some of our f{R) models. We also 
note that our analysis can be extended to scalar-tensor 
models with couplings Q of the order of 1 between dark 



and non-relativistic matter in the Einstein frame 
39(. It will be of interest to investigate the growth 
of matter perturbations and the resulting dispersion of 7 
in such theories. 
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